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Abstract 

Stochastic solutions provide new rigorous results for nonlinear PDE's 
and, through its local non-grid nature, are a natural tool for parallel 
computation. There are two different approaches for the construction 
of stochastic solutions: MacKean's and superprocesses. However, when 
restricted to measures, superprocesses can only be used to generate so- 
lutions for a limited class of nonlinear PDE's. A new class of superpro- 
cesses, namely superprocesses on signed measures and on distributions, is 
proposed to extend the stochastic solution approach to a wider class of 
PDE's. 



1 Introduction 

A stochastic solution of a linear or nonlinear partial differential equation is a 
stochastic process which, when started from a particular point x in the domain 
generates after time t a boundary measure which, when integrated over the 
initial condition at t = 0, provides the solution at the point x and time t. For 
example for the heat equation 

1 d 2 

d t u(t,x) = -g-^u(t,x) with u(0,x) = f(x) (1) 

the stochastic process is Brownian motion and the solution is 

u(t,x)=E x f{X t ) (2) 
E x meaning the expectation value, starting from x, of the process 

dX t = dB t (3) 
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The domain here is K x [0, t) and the expectation value in © is indeed the inner 
product /) of the initial condition / with the measure /i t generated by the 
Brownian motion at the t— boundary. The usual integral solution, 

with the heat kernel, has exactly the same interpretation. Of course, an impor- 
tant condition for the stochastic process (Brownian motion in this case) to be 
considered the solution of the equation is the fact that the same process works 
for any initial condition. This should be contrasted with stochastic processes 
constructed from particular solutions. 

That the solutions of linear elliptic and parabolic equations, both with 
Cauchy and Dirichlet boundary conditions, have a probabilistic interpretation is 
a classical result and a standard tool in potential theory [I] [2] [3J. In contrast 
with the linear problems, explicit solutions in terms of elementary functions 
or integrals for nonlinear partial differential equations are only known in very 
particular cases. Therefore the construction of solutions through stochastic 
processes, for nonlinear equations, has become an active field in recent years. 
The first stochastic solution for a nonlinear pde was constructed by McKean 
[3] for the KPP equation. Later on, the exit measures provided by diffusion 
plus branching processes [5] [5] as well as the stochastic representations recently 
constructed for the Navier-Stokes [2] [9] [10] [11], the Vlasov-Poisson [12] 
[T3] [15] , the Euler [14] and a fractional version of the KPP equation [16] define 
solution-independent processes for which the mean values of some functionals 
are solutions to these equations. Therefore, they are exact stochastic solutions. 

In the stochastic solutions one deals with a process that starts from the 
point where the solution is to be found, a functional being then computed 
on the boundary or in some cases along the whole sample path. In addition 
to providing new exact results, the stochastic solutions are also a promising 
tool for numerical implementation, in particular for parallel computing using 
the recently develop probabilistic domain decomposition method [T7] [TB] [T§] . 
This method decomposes the space in subdomains and then uses in each one 
a deterministic algorithm with Dirichlet boundary conditions, the values on 
the boundaries being determined by a stochastic algorithm. This minimizes 
the time-consuming communication problem between domains and allows for 
extraordinary improvements in computer time. 

There are basically two methods to construct stochastic solutions. The first 
method, which will be called the McKean method, is essentially a probabilis- 
tic interpretation of the Picard series. The differential equation is written as 
an integral equation which is rearranged in a such a way that the coefficients 
of the successive terms in the Picard iteration obey a normalization condition. 
The Picard iteration is then interpreted as an evolution and branching process, 
the stochastic solution being equivalent to importance sampling of the normal- 
ized Picard series. The second method constructs the boundary measures of a 
measure- valued stochastic process (a superprocess) and obtains the solution of 
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the differential equation by a scaling procedure. 

For a detailed comparison of the two methods refer to [20]. Here after a 
short review of the construction of superprocesses one illustrates the basic reason 
why this construction, when restricted to measure-valued superprocesses, can 
only be applied to a limited class of nonlinear partial differential equations. A 
wider class of superprocesses, namely superprocesses on signed measures and on 
distributions might provide useful stochastic solutions for some other nonlinear 
PDE's. 

2 The construction of superprocesses 

In the past, superprocesses have been constructed in the space M + (E) of fi- 
nite measures on a measurable space (E,B). Because this setting somehow 
restricts the range of nonlinear equations for which solutions may be obtained 
by superprocesses, it is convenient to use a wider framework. 

Let S be the Schwartz space of functions of rapid decrease and U C S those 
functions in S that may be extended into the complex plane as entire functions 
of rapid decrease on strips. W , the dual of U, is Silva's space of tempered 
ultradistributions [21] [22], which can also be characterized as the space of all 
Fourier transforms of distributions of exponential type }23j . Furthermore, for 
reasons to be clear later on, it is convenient to restrict oneself to the space W 
of tempered ultradistributions of compact support [22] . 

Denote now by (X t ,Po, v ) a branching stochastic process with values in U' 
and transition probability P 0m starting from time and v e Uq. The process is 
said to satisfy the branching property if given v — V\ + 

Po,u = Po,ut * Po,»2 ( 5 ) 

that is, after the branching [Xl,PQ tV1 ) and (Xf,P ^ 2 ) are independent and 
X] + X? has the same law as (X t , Po,v)- In terms of the transition operator V t 
operating on functions on U this is 

(V t f, vi + i/a) = (V t f, vx) + (V t f, v 2 ) (6) 

where e'^M = P 0tU e~ {f ' Xt) or 

(V t f,u) =-logP ,*e-<'' x *> (7) 

Underlying the usual construction of superprocesses, in the form that is 
useful for the representation of solutions of PDE's, there is a stochastic process 
with paths that start from a particular point in E, then propagate and branch 
but the paths preserve the same nature after the branching. In terms of measures 
it means that one starts from an initial 8 X which at branching originate other 8's 
with at most some scaling factors. It is to preserve this pointwise interpretation 
that, in this larger setting, we are restricting to ultradistributions in Uq. Any 
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ultradistribution in Uq may be represented as a multipole expansion at any point 
of its support, that is as a series of S's and their derivatives [22]. Therefore any 
arbitrary transition in the process X t in U' may be associated to a branching 
of paths in E and along these new paths new distributions with point support 
propagate. As a result the construction now proceeds as in the measure-valued 
case. 

In M = [0, oo) x E consider a set Q C M and the associated exit process 
£ = (£, t ,Ho, x ) with parameter k defining the lifetime. The process stars from 
x € E carrying along an ultradistribution in Uq with support on the path. At 
each branching point of the £ t — process there is a transition ruled by the P 
probability in Uq leading to one or more elements in Uq. These Uq elements 
are then carried along by the new paths of the £ t — process. The whole process 
stops at the boundary dQ, finally defining a exit process (Xq,Pq :U ) on Uq. If 
the initial v is 5 X one writes 

u{x) = (V Q f,u) = -logP ,xe- (/ ' x «> (8) 

(/, Xq) being computed on the (space-time) boundary with the exit ultradis- 
tribution generated by the process. 

The connection with nonlinear PDE's is established by defining the whole 
process to be a (£,ip) — superprocess if u (x) satisfies the equation 

u + Gq^j (u) = K Q f (9) 

where Gq is the Green operator, 

G Q f(r,x)=U 0>x f f(a,t,)da (10) 
Jo 

and Kq the Poisson operator 

K Q f(x)=Ilo, x l T<00 f(U (11) 

ip (u) means ip (0, x; u (0, x)) and r is the first exit time from Q. 

The superprocess is constructed as follows: Let tp(s, x;z) be the branching 
function at time s and point x. Then for e - ^ ' 2 ^ = Po,xe~^' XQ ^ one has 



Po, x e-V' x ^ = e- w ^ = no,, 

L v 

(12) 

t is the first exit time from Q and / (r, £ T ) = (/, Xq) is computed with the exit 
boundary ultradistribution. For measure-valued superprocesses the branching 
function would be 

oo 

¥>(s,y,z) = c^2p n (s,y)z n (13) 
o 

with YlnPn = 1 an d c denoting the branching intensity, but now it may be a 
much more general function. 
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For the interpretation of the superprocesses as generating solutions of PDE's, 
an essential role is played by a transformation of Eq.(fT2j) that uses Jl ke~ ks ds = 
1 — e~ fer and the Markov property IIo,a;l s < T n Si £ s = IIq )X 1 s<t . This is lemma 
1.2 in ch.4 of Ref.[5]. Because it only depends on the Markov properties of the 
(£ t , Ho tX ) process it also holds in this more general context. A proof is included 
in the Appendix with the notations used in this paper. 

Using the lemma, Eq. dT^I) for e~ w ( 0,x ' is converted into 



e -»(o,») = Uox 


e -f(T,t T ) 


f ft f ds 








Jo 





(14) 

Eq.© is now obtained by a limiting process. Let in (fl"4"|) replace w(0,x) by 
[3w/3 (0,x) and / by fif. /3 is interpreted as the mass of the particles and when 
the Wg-valued process Xq — > (5Xq then P M — » Pm . 



-*»(o,«) = n 



ds 



(15) 



Two scaling limits will be used in this paper. The first one, which is the one 
used in the past for superprocesses on measures, defines 

.,(!) 



1 



f (i) 

if} 



i 



-pi 



)/f> 



and 



0,x;l 



#(0,,;^)=|( V ( 
one obtains from (I15[) 

uf (0, x) + n . x f T dsipf (s, a*, 1 *) - (r, e T ) 



that is 



X /3 



G, 



(i) 





(16) 
(17) 

(18) 
(19) 



When /3 — > 0, /i 1 ' — » / and if ^ goes to a well defined limit ip then up tends to 
a limit u solution of © associated to a superprocess. Also one sees from (fl"6l 
that in the /3 — ► limit 



logPo^e' 



(20) 



as in Eq.©. The superprocess corresponds to a cloud of particles for which 
both the mass and the lifetime tend to zero. 



3 Measure- valued superprocesses and nonlinear 
PDE's 

Here one restricts oneself to measure-valued superprocesses, that is, in terms 
of paths, to S's propagating along the paths of the (£ t ,IIo,x) process and sim- 
ply branching to new S measures at each branching point. Let us construct a 
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superprocess providing a solution to the equation 



du Id 



a 



dt 2 dx 2 

for 1 < a < 2. Comparing with (|9]) one should have 

tp (0, x; u) = u a 
Then from (JTTJ) and (fTB")) . with z = 1 — /3u^ one has 



(21) 



n , /? (1)« , (! _ z ) 

ip(0,x;z) = } p n z =z + —uy =z + 



z + - J-^ f 1 - az + " (Q ~ 1} , 2 - ° (a ' 1} , (Q ~ 2) , 3 

fc^/T^ 1 V 2 3! 



(22) 



Choosing fc^ = g °_ t the terms in 2 cancel and for 1 < a < 2 the coefficients 
of all the remaining z powers are positive and may be interpreted as branching 
probabilities. It would not be so for a > 2. Then 

Po = -\ Pi=0; ••• p n = t*L{ < *} n >2 (23) 
a a \ n J 

with J2 n Pn = With this choice of probabilities p n for branching into new 
5 measures and with hp = -os=r and (3 — > one obtains a superprocess which, 
through provides a solution to the Eq. (j2"Tj) . a = 2 is an upper bound for 
this representation, because for a > 2 some of the p' n s would be negative and 
would not be interpretable as branching probabilities. 
For the particular case 

du 1 d 2 i 



.2 



U 



dt 2dx 2 " (24) 

1 2 

pi = 0; po = p 2 = -; fcp = - (25) 

When /3 — > 0, the solutions are given by ((5J) and the superprocesses correspond 
to the scaling limit of processes where both the mass and the lifetime of the 
particles tends to zero and at each bifurcation point one has probability po of 
dying without offspring or creating n new S measures with probabilities p n . 

Superprocesses are usually associated with nonlinear PDE's in the scaling 
limit (3 — > of (fTTT) - (jT5|) . However other limits may also be useful. For example 
with with p n = S n _2, /3 = 1 and hp = 1 one obtains 

^(0,x;u«) = &( v (o,x;l-/^>)-l + 



k 




A- 



a «? )3 -)9««) 



/3 

?i 2 - u (26) 
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Therefore, in this case, one is led to the KPP equation 



du 1 d 2 u , , . 

Tt = 2d^- U+U (27) 

However in this case, because /3 = 1 instead of j3 — > 0, the solution is given by 
(1 — e~ w ) instead of ([8]). Because of the natural stochastic clock provided by the 
linear u term, a stochastic solution for the Cauchy problem of the KPP equation 
may be constructed by other method [4j. However, the interpretation as an 
exit measure, allows for the construction of solutions with arbitrary boundary 
conditions. 



4 Superprocesses on signed measures and ultra- 
distributions 

Although the scaling limit (3 — > of measure-valued superprocesses allows the 
construction of solutions for equations which do not possess a natural Poisson 
clock, it has the severe limitation of requiring a polynomial branching function 
tp(s,x;z). This automatically restricts the nonlinear terms in the pde's to be 
powers of u. In addition, these terms must be such that all coefficients in the z n 
expansion in Eq.(|13p be positive to be interpretable as branching probabilities. 
As seen before, it was this requirement that led to the restriction 1 < a < 2 in 

The variable z that appears in <pg(s,x;z) is in fact z = e ^ wt ^ T s ^s> = 
Pq,xG~^^ X ' ■ When restricting the superprocess to measures, the delta measure, 
at each branching point, may at most branch into other deltas (with positive co- 
efficients) and therefore tp (s, x; z) must be a sum of monomials in z. When one 
generalizes to UL\ ultradistributions of point support^, changes of sign and tran- 
sitions from deltas to their derivatives are allowed. In the end, the exponential 
e —{Pf,x) w jv| kg computed by evaluation of the function on the ultradistribution 
that reaches the boundary. To find out the equation that is represented by the 
process one needs to compute the ip g (0, x; up) of Eq. (fT7|) for the corresponding 
if (s, x; z) in the /3 — > limit. Recalling that tp (s, x; z) — tpp (s, £ s ; e -,9 ™'(' r ~ s '^') 
and z = e~P Wl3 , one concludes that there are basically two new transitions at 
the branching points: 

1) A change of sign in the point support ultradistribution 

e </3/,<5*> = e /3f{x) _^ e <J3f,-S a ) = e -fSf{x) £ 2g ) 

which corresponds to 

z -> - (29) 

and 

1 Because distributions of point support are a finite sum of deltas and their derivatives [23], 
one could have considered only distributions of point support rather than compact support 
ultradistributions ML. However in ML one is not restricted to finite sums. 
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2) A change from 5 (n) to ±5 {n+1) , for example 

e WW = -> e< w > ±4 -) = e^'^ (30) 

which corresponds to 

Case 1) corresponds to an extension of superprocesses on measures to superpro- 
cesses on signed measures and the second to superprocesses in Uq. 

How these transformations provide stochastic representations of solutions 
for other classes of pde's, will be illustrated by two examples: 

First, let 

ip (1) (0, x; z) = Pl e d * log z + p 2 e- 9 * log z + p 3 z 2 (32) 

This branching function means that at the branching point, with probability p\ 
a derivative is added to the propagating ultradistribution, with probability pi a 
derivative is added plus a change of sign and with probability p 3 the ultradistri- 
bution branches into two identical ones. Using the transformation and scaling 
limit (I16p one has, for small f3 



(33) 

z -> z 2 = ( 1 - Pu { g } ) = 1 - 2(3u { g 1] + /3 2 u { f p 2 (34) 



3\ 



Then, computing ip g ^0, x; u^^j with Pi — P2 — j and £>3 = \ one obtains 
tf«(o, «; uf) = |(^)(0,x;z)-z) 

(o,,;i-/^)-i + /3u«) 

-f{d x uf)\ l -Puf* + 0{?)) (35) 



/3 

fc/3 /l 



/3 V 8 

meaning that, with kp = 4, the superprocess provides, in the /3 — >■ limit, a 
solution to the equation 

du ld 2 u 2 1 2 / Q <A 

W = 29^- 2M -2 ( ^ M) (36) 
For the second example a different scaling limit will be used, namely 



uf = ^{e^-e-^) ■ ff> = ±(e?f-e-ef) (37) 

,( 2 ) _^ .... ar .A f( 2 ) 



Notice that, as before, u^' —¥ wp and fp—>f when j3 —> 0. In this case with 
z = e^ Wf> one has 



-2$uf + 2^Jj3 2 uf )2 + l 

2 - 2/3uf + p 2 uf 2 + O (/3 4 ) (38) 
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and 

\ = 2?uf + 2^ 2 uf 2 + 1 

= 2 + 2(3uf + (3 2 uf 2 + O (/3 4 ) (39) 
For the integral equation, instead of (|T8|) . one has 

uf (0, x) + H , x jT ds4 2) (s, £ s ; uf) = U^ x ff (r, £ T ) (40) 

with 

V^ 2) (0, x- uf) = k p CL L (0, x;z)-<p U x; ±\ \ - ufA (41) 

Let now 

V {2) (0,x;z) =piz 2 +p 2 - z (42) 

This branching function means that with probability p\ the ultradistribution 
branches into two identical ones and with probability p 2 it changes its sign. 
Therefore, in this case, one is simply extending the superprocess construction 

to signed measures. Using ([38} and (|39]) one computes ijjf ^0, x; uf) obtaining 

4 2) (0, uf) = k { - Pl 8uf (l + \f?uf) + V.uf - uf + O } 

. (43) 

and with Pi = j^',P2 = jq and hp = one obtains in the in the j3 — > limit 

4 2) (0,x;4 2) ) ->-uf 3 (44) 

meaning that this superprocess provides a solution to the equation 

du 1 <9 2 u 



dt 2dx 



u 



(45) 



In conclusion: Extending the superprocess construction to signed measures 
and ultradistributions, stochastic solutions are obtained for a much larger class 
of partial differential equations. 

5 Appendix: Proof of the lemma 

Let 

u(x,t) = n o ,Je- kt u(a t ,0)+ [ ke- ks $(l; s ,t-s)d S } (46) 



(i 
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Then 



U °,*J q ku{£ s ,t-s)ds = n Q , x jjf ke- k ^u^ s+t _ s ,0)ds 

+ kds kds'e- ks '<S>(£ s+s ,,t- s- s')| 



(47) 



Summing @5) and 07]) 

u (x, i) + Hq,x / ku (£ s ,t — s) ds 
Jo 

= U^ x ^e- kt +^ ke-^-^ds^jui^O) 

+k [ e- ks $(£ s ,t~ s)ds + k f ds f kds'e~ ks '^^ s+s ,,t-s-s')ds'\ 
Jo Jo Jo ) 



(48) 



Changing variables in the last integral in (|48p from (s, s') to (s, cr = s + s') one 
obtains for the last term 



and finally 



k I da I kdse-~ k ^~ s ^{£, a ,t-a)ds 
Jo Jo 

u (x, t) + Uo.xk / u(£ s ,t — s)ds 
Jo 

= IIa x j u (£ f ,0) + fc^ $(£ s ,t-s)ds 



(49) 
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